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ABSTRACT 

^ I The anomalous bump in the cosmic ray positron to electron ratio at 10 GeV 

O i can be explained as being a component from a point source that was originally 

^ I harder than the primary electron background and degrades due to synchrotron 

^ ' and inverse Compton losses in the Galaxy while propagating to the Earth's vicin- 

ity. The fit is better than can be obtained with homogeneous injection and is 
attributed to a minimum age threshold. Annihilating neutralinos can provide 
a fair fit to the data if they have a mass just above 1/2 the mass of the Z° 
\Q • and if they annihilate primarily in distant density concentrations in the Galaxy. 
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A possible observational consequence of this scenario would be intense inverse 
Comptonization of starlight at the Galactic center, with a sharp energy cutoff in 
the emergent photons as a possible signature of the neutralino mass. 

Subject headings: cosmology: dark matter - diffusion - elementary particles - 
galaxy: center 



1. Introduction 

The possibility that weakly interacting dark matter particles (WlMP's) could annihilate 
into detectable cosmic radiation was suggested by Silk & Srednicki (1984). Tylka & Eichler 
(1987) noted a reported positron excess, curiously localized near 10 GeV [Mueller & Tang 
(1985, 1987); Barwick et al. (1997, 1998)] and considered whether it could be due to the 
annihilation of photinos (as a simple example of neutralinos) in the tens of GeV mass range. 
The difficulty was that this process, given the laboratory constraints on the neutralinos, 
seemed to fall short of providing enough positrons, and the results were not published. 
Various papers on this excess eventually appeared [Tylka (1989); Eichler (1989); Turner & 
Wilczek (1990); Coutu et al. (1999)], and some noted that the potential for positron excess 
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could be bolstered by dumpiness in the annihilating dark matter or by decay of weakly 
unstable dark matter particles. 

The approach usually found in present literature is to try and fit the overall e+/ (e+ + e~) 
ratio, without giving special attention to the curious behaviour at 10 GeV, see Hooper & 
Silk (2004). Baltz & Edsjo (2001) considered a whole class of minimal standard supersym- 
metric models and failed to get any non-monotonicity in the e~^/{e'^ + e~) ratio. Eichler & 
Maor (2004) (henceforth Paper I) considered annihilation of particles through the Z°-channel 
(which we shall henceforth refer to as virtual Z°-decay) noted that non-relativistic virtual 
Z° decay (i.e. when the rest mass of the annihilating dark matter particle is slightly above 
1/2 the Z" mass) provide a remarkably good fit to the observed e+/(e^ + e~) ratio below 
10 GeV mainly due to the positrons that emerge from decaying muons. At higher energies, 
however, the predicted e'^/(e'^ + e'~) ratio rose above the observed values within conventional 
assumptions about the injection and propagation. In particular, it was assumed in Paper I 
that the positrons and primary electrons are each injected with the same spatial profile, and 
that their propagation in the Galaxy is identical. The reason for this rise is that some Z°^s 
decay directly into high energy e~^e~ pairs so that the e"*" energy is half the Z° mass, and 
this gives rise to a high energy bump in the e+/ (e+ + e^) ratio at about 50 GeV. While this 
bump can be partially washed out by losses and escape, it was found that the high energy 
e"''/(e+ + e") ratio is nevertheless apparently too high to fit the observations to within la- 
error bars. As discussed in Eichler (1989) this is a generic problem for any positron source 
that is significantly harder than the primary electrons above 10 GeV. 

However, the dark matter annihilation scenario for explaining the positron excess in 
any case requires that clumping of the dark matter, and a likely place for this is near the 
Galactic center. This means that the positrons in our neighborhood that are dark matter 
annihilation products would have a minimum age, i.e. the time needed to diffuse from the 
source to our neighborhood, and the age distribution of the positrons that make it to the 
Earth's vicinity contain fewer young positrons than the age distribution that one associates 
with the standard leaky box model. In this letter we consider that the positrons are injected 
by an effectively point source at a finite distance, and show that it greatly improves the 
fit over that obtained in Paper I. The positron bump at ~ 10 GeV can be attributed to 
halo-type age ~ 3 x 10^ yr for the positrons, for over such a lifetime, positrons losing energy 
by synchrotron and inverse Comptonization would end up at about this energy. 

We will find that obtaining a good fit from a single source with a single diffusion coeffi- 
cient is difficult. However, it is well known that below a certain intensity, e.g. far enough up 
front in a diffusion front, cosmic rays can freely stream. Such behavior is observed upstream 
of the Earth's bow shock. Theoretical reasons for such free streaming include the difficulty 
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of resonantly scattering cosmic rays through 90 degree pitch angle in the linear wave am- 
plitude regime. Also, a larger, counterstreaming , finitely stable component of cosmic rays 
would stabilize the smaller free streaming positron component. We therefore allow for the 
possibility that a small fraction of the positrons freely stream, and arrive at the Earth's 
vicinity much younger than the rest. We find that this improves the fit still further. We 
find that the low energy non-monotonicity, which appears in the injection spectrum from 
non-relativistic Z° decay Paper I, can be produced also by using a simple power law for 
the injected spectrum (using the same propagation model). It can thus be produced by the 
combination of a harder (but monotonic) spectrum of injected positrons and propagation 
effects. We conclude that it is still too early to unambiguously interpret the low energy 
behavior of the spectrum as a signature of self annihilating dark matter. 



2. Equations and Results 

The steady state diffused equation for the particle number density, n{x,r), is 
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X — E/mz, -D is a diffusion operator, R — Bx is the escape rate with B ~ few x 10 



= -^x^ with A — 8.5 x 10~^^ erg/s is the Compton loss rate, corresponding to an 
electromagnetic energy density in the Galaxy of 10~^^ erg/cw?, L is a distance scale, and 
I{x) is the spectrum injected by a point source. 

We assume as in Paper I that the primary electrons and background positrons are 
injected homogeneously, Dri}, — 0: 
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For the Z" decay injected spectrum we take the diffusion to be one-dimensional with 
a diffusion coefficient V, D — 'D-^, and with boundary conditions such that ^\r=L — 
(conserving the number of particles except for the escape term). With these boundary 
conditions, the solution to eq. (1) is 
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While V (the diffusion coefficient) and L (the size of the leaky box) are free parameters, we 
took r = 8 Kpc, the distance to the galactic centre. K = V/r^ gives the inverse time for 
diffusion. 

We have chosen a one dimensional diffusion because it gives somewhat better results than 
3 dimensional diffusion. This is physically plausible if one considers magnetic fields which 
will confine the movement of the charged particles. So the geometry is tube-like, with an 
effective cross section such that the total volume is the galactic volume, (20 kpc)^. 

Iz{x) is the Z decay products, 

Iz{x) ^ N 0.0344/e(x) + 0.0344/^(x) + 0.0069/^(x) + 0.69164(x) ) 

1 

2 



(5) 



Ie(x) 


= M 


x 




2 
3 




Ir{x) 


2 
3 





= 3(5- 2>Qx^ + 2>2x^ ) 
?6x^ + 32x^ ) 



+ 



95 



108x^ + 

3 



h{x) 



dx 



2 
9 

14 

"9 _ 

ai — 3, bi — 10 
02 2, 62 = 4 



1408 



x^ - (25 + 324^2 + 128x^) In {2x) 



f 



-10 



X 



< X < 0.1 
O.K X < 1 



Each Ich describes the ch channel of decay, and the pre-factors correspond to the branch- 
ing ratios. The calculation was done in zeroth order, assuming 3 massless families and 
neglecting the top quark, for details see Paper I. Following the discussion there, we take 
= 1.3 X 10~^^ l/(cm^ s) as the annihilation rate per unit volume. 



Fig. (1) shows a fit with a single point source of Z° decay and a single diffusion 
coefficient. The good fit to the low energies from Paper I is still present, but at the price 
that the excess in energies toward x = 1/2 is now is suppressed by the finite age effect. As 
the figure shows, we are now facing a scenario which is opposite to Paper I; the finite age 
effect tends to suppress the high energy excess at the price of kiUing it off altogether. 
However, there are several possibilities that avoid this problem: There may be more than 
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one source, and there may be more than one route (roughly guided by magnetic field lines) 
by which the particles diffuse or freely stream from the source to our vicinity. High energy 
particles diffuse much less than low energy ones because they are fewer in number and create 
less waves. So their self-generated scattering is less efficient. Thus, the fraction of free 
streaming particles should be higher at higher energy. Fig. (2) shows a combination of two 
Z° decay components, an older, larger one that arrives via diffusion, and a younger, smaller 
component that has managed more free streaming. This figure illustrates that if one takes an 
age distribution into account, the fiexibility in adjusting the high energy spectrum becomes 
much larger, and can be fitted to the data. 

For sake of comparison, we also include a power law injected spectrum, fig. (3) shows 
various power laws, and fig. (4) shows two components with different ages. We find that as 
long as the injected power law is hard enough, one can produce a low energy (5 — 10 GeV^ dip. 
The quahty of the fit is almost as good for a power law as for virtual Z° decay. We consider 
the low energy dip to have qualitatively more significance than the higher points and have 
emphasized those data points accordingly in choosing the best fit. We have deliberately not 
quantified this with the standard statistical measures. Trying to get the statistically best 
parameters (for either power law or Z° decay as injected spectrums) would wash out the low 
energy behavior that we are focusing on. 

Although we can reproduce the 7 GeV dip, the peak at ~ 15 GeV is still too big for 
the HEAT data (though too small for the earlier data). This seems to be a generic feature 
of our results, regardless of whether the injection source is virtual Z° decay or a power law. 
The problem would be worse if the virtual Z° had an energy well above mz- 

3. Possible Observational Consequences 

The hypothesis that the neutralino mass is only slightly more than half the Z° mass 
is motivated by several factors: The annihilation cross section can be resonantly enhanced 
by a factor of 2 or 3 more than that during annihilation in the early universe, when the ther- 
mal broadening of the Z resonance somewhat exceeded its natural width, Greist and Seckel 
(1991). Moreover, assuming the smallest allowable mass allows the greatest annihilation rate 
since the annihilation reaction rate is fixed by the condition that it allows a given cosmic 
dark matter contribution. (Although dark matter clumping can enhance the annihilation 
rate, a plausible level of such enhancement is limited by observational constraints on dark 
matter clumping that are set by stellar distributions in galactic centers.) Making just 
above mz causes the Z° resonance to be asymmetric, but this would be hard to measure 
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experimentally because of the weak coupling of the emerging neutrahnos at CM coUision 
energies above ^rriy.. On the other hand, that the annihilation cross the virtual Z° must be 
close to its mass shell if it is to provide a decent fit suggests that its loop corrections would 
be large and it might be discernable or falsifiable with particle collider data on processes 
that depend on such loop corrections. 

In an astronomical context, a possible observational consequence of a point source of 
positrons at the Galactic center could be inverse Comptonization of starlight, which is far 
more intense than at a typical point in the Galaxy. The profile of Galactic starlight near 
the Galactic center is given by Kent (1992). The derived photon energy density is then 
U{r) — 4.3 X 10~^(r/pc)~°'^^ erg/s. Assuming the positrons are produced within a typical 
radius r of the Galactic Center, they produce a minimum of -B/cCt) = J^^^'^ ^^o-TU{r)dr 
in inverse Compton (IC) scattered starlight before escaping the central 100 pc region, and 
the luminosity L{> '-jrrif.c^) above 7^6^/1, where is the typical energy of the pre-scattcred 
starlight photons, is I{j'meC^)Eic{'~i')d'~^' . The most energetic e^e^ pairs alone, that is 
those which result directly from the Z° decay {Ie{x) in eq. 6), will produce IC luminosity of 
3.8 X 10^5 erg/s. 

Fig. (5) shows the logarithmic derivative of the IC luminosity due to the positrons 
only, — 7^L^2 = — (iL/(iln(7^), as a function of the square of the positron Lorentz factor, 
7^. Shown in the figure is the minimum IC luminosity as a function of the frequency scaled 
to the frequency of the prc-scattcrcd photon (x axis). The minimum luminosity assumes 
that the positrons emerge from the central region in a straight line. If the mean free path 
A is less than 100 pc, then the predicted luminosity goes up by roughly a factor of (lOOpc/A). 

For dark matter annihilation that yields direct monochromatic e+e^ pairs at Lorentz fac- 
tor 7o, this would translate into a sharp cutoff in the IC gamma rays of 7oep/i or about 37^ eV . 
In our particular example m-^ ~ rnz/'i, this would lead to a cutoff at 10^°ep/i ~ 30 GeV. This 
would in principle be detectable by MAGIC, see for example Cortina (2004), if the location 
were suitable for observing the Galactic center. Alternatively, it could be detected by HESS 
if the energy threshold could be pushed to below 30 GeV. This scenario would not explain 
the TeV photons from the Galactic center recently reported by the HESS collaboration, 
Aharonian et al. (2004). If, however, there is annihilation in the Galactic Center of heavier 
dark matter particles, then direct e^e~ pairs might be detectable via such a cutoff in the 
TeV gamma ray spectrum at ra^jl. 

In conclusion, we find that the cosmic ray positron data can be fit with more than one 
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hard source of positrons provided that a) they have a chance to lose energy before escaping 
the Galaxy and b) they have a minimum age (e.g. they come from discrete, distant sources), 
unlike the background primary electrons. They need not be from dark matter annihilation, 
but a best case scenario for this is not confidently ruled out by existing data. Detection of 
inverse Compton radiation with good energy resolution can in principle provide information 
as to the spectrum of the positrons and their point of origin. 
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Fig. 1. — The /{e^ + e ) as a function of a; = E/mz, for a single source Z decay injected 
spectrum. A = 8.5x10-^^ r = 7.1x10-^^ -,C = 4.0x10-^9 -J^ D = 1.3x10'^^ 



and X = 1.9 X 10"^^ \. Data taken from Barwick et al. (1997) (black) and Mueller & Tang 
(1987) (grey). 
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Fig. 2. — The e+/(e+ + e ) as a function ol x — E/mz, for a combination of 2 sources of 
Z decay injected spectrum. ^ = 8.5 x 10"^^ ^, B ^ 7.6 x 10"^^ -, C = 4.9 x IQ-^^^ 
D = 1.3 X 10-^^ Ki = 2.8 X 10-^^ ^ and K2 = 2.8 x 10"^^ i. The ratio between the 

two components is 1 : 5. Data taken from Barwick et al. (1997) (black) and Mueller & 
Tang (1987) (grey). 




Fig. 3. — The e"''/(e+ + e ) as a function oi x = E/ruz, for various power laws, Nx"^, as 
the injected spectrum. A = 8.5 x 10"^^ B = 4.4 x 10"^^ C = 1.7 x IQ-^^ 
D = 1.1 X 10-^^ and X = 6.6 x 10"^^ Data taken from Barwick et aL (1997) 

(black) and Mueller & Tang (1987) (grey). 
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Fig. 4. — The /{e^ + e ) as a function of x = E/mz, for a combination of 2 sources 
of power law {w = -0.3) injected spectrum. A = 8.5 x 10"^^ ^, S = 4.4 x 10"^^ ^, 
C = 1.7 X 10-2^^ D = 1.1 X 10-=^^ = 7.2 X 10-^^ i and K2 = 1.2 x IQ-^^ K 

cm'^ s ' cm" s' ^ s s 

The ratio between the two components is 10 : 1. Data taken from Barwick et al. (1997) 
(black) and Mueller & Tang (1987) (grey). 
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Fig. 5. — Differential IC luminosity due to the positrons only, —dL/d\n{'-f'^) x (20 kpc)^/V 
as a function of 7^. 



